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In this Comment, we show that the experimental data reported by Li et al. (Phys. Rev. B 
77, 134501 (2008)) do not support the phonon specular reflection at low temperatures for high-Tc 
cuprates, because the phonon mean free path £ calculated from their thermal conductivity data is 
much smaller than the averaged sample width above ~ 100 mK for both Nd2Cu04 and YBa2Cu306.o 
single crystals. 

PACS numbers: 72.15.Eb, 74.25. Fy, 74.72.-h 



In a recent paper, Li et al^ reported their 
low-temperature thermal conductivity (k) studies on 
Nd2Cu04 single crystals. Based on their examinations 
of the effect of sample size and surface roughness on the 
phonon heat conductivity Kp, they concluded the domi- 
nance of specular phonon reflection in high- Tc cuprates 
and proposed an empirical power law T" (a < 2) to de- 
scribe Kp/T. They further used this empirical formula to 
re-analyze our published data^ for Zn-doped YBa2Cu30y 
(YBCO) and claimed that the universal thermal conduc- 
tivity of d-wave quasiparticles at T — > was confirmed. 
In this comment, we show that both their experiments 
and data analysis are problematic or self-inconsistent, 
and therefore their conclusion on our data of Zn-doped 
YBCO lacks reasonable ground. 

Li et al.^ showed that the thermal conductivities of 
a wide Nd2Cu04 sample and a narrow one (cut from 
the wide sample) differ, as did the thermal conductiv- 
ities of an as-grown sample and the same sample after 
surface roughening, and claimed that the observed dif- 
ferences support their assertion that phonons are in the 
boundary scattering regime below 0.5 K and that the 
specular reflection off sample surface is important, caus- 
ing the (a < 2) dependence of k in a rather broad 
temperature range. Let us first point out that, contrary 
to what Li et al. claimed, the phonons are not really 
specularly reflected in their samples by correctly calcu- 
lating the phonon mean free path £ from their data and 
comparing it with the averaged sample width W , which 
is usually taken to be 2/-\/7r times the geometrical mean 
width w.^ Note that previous studies on the specular re- 
flection unanimously found'^''^'^ that £ gets much longer 
than W when the specular reflection becomes important. 
Actually, since Kp is equal to ^Cvi where C = (3T^ is 
the phonon specific heat and v is the averaged sound ve- 
locity, and both /? and v are known experimentally for 
Nd2Cu04,— i^i^ one can easily calculate £ in Li et a/.'s 
samplesii The coefficient /3 of Nd2Cu04 were reported 
to be about 0.42-0.50 mJ/mol-K.^*i^ To avoid underesti- 
mating the value oi £, we choose (3 — 0.42 mJ/mol-K for 
the calculation. Another parameter, the averaged sound 
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FIG. 1: (color online) Temperature dependence of the calcu- 
lated phonon mean free path £ divided by the averaged sample 
width W for Nd2Cu04. The raw data of n are taken from Ref. 



velocity, can be calculated hy v = Xli Si 
where i represents the three acoustic phonon modes. 
Here, we use the experimental values of the in-plane 
sound velocities of Nd2Cu04 from Ref. ^91, vl = 6050 
m/s, vti — 4220 m/s and vt2 — 2460 m/s (which are 
also used by Li et al.^), to obtain v = 2940 m/s. Figure 
1 shows the temperature dependence of £ calculated from 
the data in Ref. Ij. One can see that for most cases £ only 
becomes comparable to W at the lowest temperature and 
it stays much smaller than W in the temperature range 
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where Li et al. argued^ that the boundary scattering is 
taking place (below 0.5 K). For surface-roughened sam- 
ple, £ is so small that it is much smaller than W even at 
the lowest temperature. Therefore, the magnitude of £ 
is obviously inconsistent with the assertion that phonons 
are boundary scattered below 0.5 K in Nd2Cu04. In fact, 
since Kp = ^{/3T^)v£, a departure of Kp from the de- 
pendence always comes from a T dependence of and 
it occurs either when the specular reflection is impor- 
tant (which causes i > W) or when phonons are not in 
the boundary scattering regime (in which case £ < W). 
Therefore, one can conclude that in Nd2Cu04 at tem- 
peratures higher than ~100 mK, phonons are not in the 
boundary scattering regime and the phenomenological 
T"+i (a < 2) dependence comes from some additional 
scattering that only reduces £M' 

Furthermore, one can see that a comparison of the 
thermal conductivity data between the wide sample and 
the narrow one indicates some self-inconsistency if one 
tries to employ the specular reflection to explain the 
T""*"^ dependence: When the sample is gradually cooled, 
one expects that £ gradually increases when the micro- 
scopic scattering is becoming less effective and finally £ 
is equal to the averaged width, where the phonons en- 
ter the boundary scattering limit (about 0.5 K in Li et 
aL's speculation). Upon further cooling, £ would become 
larger than the averaged width if the specular reflection is 
at work. Naturally, the characteristic temperature where 
£ is becoming comparable or equal to the averaged width 
should not strongly decrease with the sample size; rather, 
it should move to higher temperature upon decreasing 
the width. This is clearly opposite to what Li et a/.'s 
data imply, as is demonstrated in Fig. 1(a). Another self- 
inconsistency appeared when Li et al. compared the data 
between the as-grown sample and the surface-roughened 
onefi that is, the two data sets differ by as much as a 
factor of 2 at 0.5 K where k should already be mostly 
governed by microscopic scattering mechanisms, which 
should give a sample-size-independent thermal conduc- 
tivity. (Here, we are not talking about thermal con- 
ductance.) This is clearly inconsistent with what they 
claimedii Note that the cutting or sanding process could 
bring some cracks into the samples, and our own expe- 
rience taught us that small cracks can easily be created 
when flux-grown crystals are cut; the narrower and thin- 
ner the sample, the more significantly the cracks affect 
the transport results. 

It should be noted that Li et al. showed an estima- 
tion of £ for the surfaced-roughed sample by using the 
following formulafi 

np = ^7:'kB{^f<v-'>£, (1) 

in which only one parameter < > is necessary for 
the calculation. To obtain this formula, there is an ap- 
proximation for the phonon specific heat coefhcient 
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FIG. 2: (color online) Temperature dependence of the calcu- 
lated phonon mean free path £ divided by the averaged sample 
width W for YBa2Cu306.o. The raw data of k are taken from 
Ref. la. 

where < > is the average of the inverse third power of 
the long-wavelength phase velocities of the three acoustic 
modes.— 1^ Calculating /3 with this approximation using 
measured sound velocities would be valid if the sound ve- 
locity solely reflects the dispersion of the acoustic phonon 
near the F point. However, Eq.([T]) is apparently not pre- 
cise when one compares the experimental results of f3 and 
the sound velocities j^ii^ Using the experimental velocities 
from Ref. d, one obtains /3 = 0.20 mJ/mol-K from for- 
mula 12]), which is smaller than the experimental value 
of (3 (0.42-0.50 mJ/mol-K)4 In this sense, Li et al.'s cal- 
culation using formula ([!]) is not well supported by the 
experimental literature. Furthermore, their calculatioi>i 
itself is problematic and £ is at least 30% overestimated: 
they used a roughly averaged value of the velocities {v 
— 4000 m/s) for calculating < v~'^ > through l/w^ {— 
6.25 X 10~® s^/m^); however, one can easily notice that 
< w"^ > should be the averaged value of the inverse sec- 
ond power of the sound velocities of the three acoustic 
modes, which is < > — 8.29 x 10~* s^/m^. 

Therefore, Li et aZ.'s datai^ for Nd2Cu04 do not re- 
ally support their speculation that the specular phonon 
reflection is important in cuprates. Then, what about 
another cuprate, YBCO, which they simply assumed to 
have the same phonon transport properties as Nd2Cu04? 
One can easily find the answer from the thermal conduc- 
tivity data of an insulating YBCO sample with y = 6.0, 
reported by the same group. ^ In Fig. 2, we show the tem- 
perature dependence of £ in this sample where the heat 
transport is purely phononic, using the data from Ref. 
@. It is clear that the YBCO crystal shows a standard 
phonon transport behavior, with the boundary scattering 
limit only achieved below 140 mK, where £ is close to the 
averaged size and k shows a dependence. Actually, it 
was already demonstrated^ that the low-T k/T data of 
this sample below 140 mK can be reasonably well fit to 
a + bT^ with zero intercept (a = 0). 

Neglecting those circumstances, Li et al} used the for- 
mula k/T = a + bT" to re- fit our low-T thermal conduc- 
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tivity data of the superconducting YBCO and Zn-doped 
YBCO [Ref. g and concluded that kq/T is unchanged 
upon Zn doping. Based on such an analysis, they claimed 
that our data confirm the universal thermal conductivity 
of d-wave quasiparticles for T — > 0. However, it must be 
pointed out that even if the phonon heat conductivity 
showed a T"+^ dependence, the fitting to k/T = a + bT" 
is meaningful only when the electronic heat conductivity 
is a linear function of temperature. In this regard, one 
should not forget an important fact, which was first re- 
ported by the Taillefer group^^ that the electronic con- 
tribution to the heat transport in clean YBCO increases 
rapidly with temperature following a law, due to the 
thermal creation of heat-carrying quasiparticles near the 
nodes. If one takes into account this dependence of 
the electronic term at T 7^ 0, one can easily see that 
the analysis of our YBCO data by Li et al. is obvi- 
ously erroneous; in contrast, our original fittingSi with 
k/T = a + bT^ is robust against the appearance of the 
electronic contribution and, hence, is reliable for ex- 
tracting the residual thermal conductivity. In fact, the 
"main effect" of Zn doping to cause a dramatic suppres- 
sion of the slope of k/T, which Li et al. noted in Ref. 1, 



is mostly due to the diminishment of this electronic 
term. 

It would also be useful to note that Nd2Cu04 used 
in Li el aVs work is not an ideal system for study- 
ing the purely phononic heat transport because of its 
complicated magnetism In fact, theoretically the Nd 
magnon modes are predicted to be effectively gapless in 
zero field, so the conjecture of Li et al. that the Nd 
magnons are gapped and do not contribute to the heat 
transport does not have any theoretical support. If the 
magnon contribution is not negligible in Nd2Cu04, the 
true mean free path of phonons would be even smaller 
than those shown in Fig. 1. 
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